Calibration robust entanglement detection beyond Bell inequalities 
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In its vast majority entanglement verification is examined either in the complete characterized 
or totally device independent scenario. The assumptions imposed by these extreme cases are often 
either too weak or strong for real experiments. Here we investigate this detection task for the inter- 
mediate regime where partial knowledge of the measured observables is known, considering cases like 
orthogonal, sharp or only dimension bounded measurements. We show that for all these assump- 
tions it is not necessary to violate a corresponding Bell inequality in order to detect entanglement. 
We derive strong detection criteria that can be directly evaluated for experimental data and which 
are robust against large classes of calibration errors. The conditions are even capable of detecting 
bound entanglement under the sole assumption of dimension bounded measurements. 



I. INTRODUCTION 

Entanglement is the most striking phenomenon in the 
quantum viforld. It provides the resource for fascinating 
new apphcations such as quantum computing, teleporta- 
tion or unconditional secure communiction. These pos- 
sibihties drasticahy sparked the interest for this resource 
and many current experiments strive to reaUze strong 
and robust entanglement. Consequently many different 
detection methods have been developed in recent years, 
for reviews see Refs. [JQ. 

Rehable entanglement verification must fulfil certain 
criteria Q ■ Most importantly it should not depend on the 
preparation procedure and the only available information 
about the presence of entanglement should be obtained 
via measurements of the underlying system. However 
there is still one open choice left: Usually each classical 
outcome is associated with an operator describing the 
measurement apparatus, say outcome "k" corresponds to 
the measurement operator M^. Equipped with this quan- 
tum mechanical meaning one can eniploy for instance the 
tool of entanglement witnesses [J-|6[ to decide the en- 
tanglement question. This standard scenario might be 
too optimistic for applications because it crucially relies 
on the correctness of the employed operator description 
of the measuring device and clearly if the true measure- 
ment apparatus functions are different then anything can 
go wrong. For mere entanglement detection these devi- 
ations might be called systematic errors but for appli- 
cations where the presence of entanglement is essential, 
secure communication being prominent example, these 
deviations are undesired pitfalls 0, Q ■ Thus in contrast 
to the completely characterised scenario there is also the 
other extreme where one does not need any specific quan- 
tum mechanical model at all. Though surprising at first 
even in this totally pessimistic device independent case 
it is possible to infer entanglement for good enough data, 
for example using Bell inequalities p4ll|. 

Detection of entanglement in a complete device inde- 
pendent manner has recently attracted a lot of interest; 



in particular verification in the multipartite settings [l2| - 
[T^ l since it was realized that this task differs from the 
corresponding non-locality one [l^ . This was surpris- 
ing because device independent entanglement detection 
and exclusion of a local hidden variable model are equiv- 
alent problems in the bipartite case [l^, [l^. Steering 
inequalities are entanglement detection methods in a hy- 
brid scenarios, i.e., one party is complete characterised, 
the other totally uncharacterised • Only a few results 
and techniques have addressed the detection of entangle- 
ment in a partially characterized setting so far. Refer- 
ences [H, [l^ considered the case of sharp, orthogonal 
qubit measurements and showed that much more entan- 
gled state can in fact be detected than with the corre- 
sponding Bell inequality. For instance it was proven that 
the bound appearing in the famous Bell correlation term 
of the CHSH inequality [3 

can be actually reduced from 
2 to when the measurements satisfy this extra con- 
straint. These results have been extended in order to 
provide even quantitative bounds on the amount of en- 
tanglement in Ref. [20| . In order to devise more robust 
entanglement detection methods that avoid fake entan- 
glement detection under misspecification of the employed 
observables a technique called squash model is very use- 
ful dll, . Applied to entanglement verification, this 
notion, usually common in quantum key distribution, can 
even be enlarged [l^ . Finally, Ref. [23| introduces types 
of Bell inequalities that can be applied if the commutator 
of different measurement settings is known. 

The purpose of the current manuscript is to investigate 
the verification task in the intermediate regime where one 
possesses some partial knowledge of the employed mea- 
surement devices. Despite its elegance the complete de- 
vice independent setting suffers from the drawback that 
it is hard to address experimentally. In fact all Bell ex- 
periments so far were not fully device independent due 
to certain loopholes, but which are often irrelevant if 
certain other assumptions hold. However such assump- 
tions like the fair sampling condition effectively can be 
regarded as partial information, e.g., that the "inconclu- 
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sive" measurement outcome is the same for both settings 
for the fair samphng case. But clearly if these additional 
assumptions are not satisfied, the corresponding imple- 
mentation of the Bell test also does not provide a positive 
entanglement check [H, [2^. Besides, the complete de- 
vice independent scenario is often also too pessimistic. 
Of course for applications like quantum key distribution 
this very pessimistic viewpoint is legitimate but for the 
mere verification of entanglement in an experiment this 
seems like breaking a butterfly on a wheel. Although 
we do not address the question how such partial knowl- 
edge can be obtained or justified, we nevertheless believe 
that certain deviations in a measurement description are 
more stringent than others. For example the assumption 
that the measurement of the electronic state of an ion in a 
trap is very well described by a qubit measurement seems 
much easier to assure than the assumption that the per- 
formed measurements are really orthogonal or that they 
are true projectors. Independent of this discussion which 
scenario is now more reasonable for which situation, our 
investigation shed some light on the question which as- 
sumptions are more crucial than others in order to verify 
entanglement. In addition, the derived entanglement cri- 
teria are more robust against calibration errors while they 
still keep a large detection strength, in particular when 
compared to Bell inequalities. 

In the following we first focus on the scenario inves- 
tigated in Ref. [l9| and analyse entanglement verifica- 
tion in the simplest possible setting of two different di- 
chotomic measurement settings per side. We draw our 
attention to different assumptions about qubit measure- 
ments and distinguish three different classes: sharp, or- 
thogonal or the completely uncharacterised qubit mea- 
surements. We provide a solution in terms of the singu- 
lar values of a corresponding data matrix and find that 
one detects a much larger fraction of states than with 
the complete device independent setting. This already 
provides examples where one detects entanglement al- 
though the corresponding Bell inequalities, the CHSH 
inequalities [l^l in this case, are not violated with the 
sole extra assumption that the dimension of the under- 
lying quantum system is fixed. Additionally we con- 
sider specific observations where the additional knowl- 
edge of sharpness or orthogonality is irrelevant for the 
detection strength and already the dimension restriction 
suffices to verify exactly the same amount of entangled 
states as with completely characterized, i.e., sharp and 
orthogonal measurements. After considering these vari- 
ous scenarios for two qubits we focus on the extensions to 
more dichotomic measurements with completely unspec- 
ified measurements restricted only by the underlying di- 
mension. We derive a criterion that is applicable for any 
of these settings and show that it is capable of detecting 
entanglement in data originating from bound entangled 
states. Moreover the criterion even shows that with un- 
characterised qutrit measurements one can verify more 
entanglement than with corresponding Bell inequalities. 
Note that dimensions d > 4 arc not relevant, because 



then one effectively equals the detection strength of the 
Bell inequalities [l^ . 

The outline is as follows: In Sec.|ll]we precisely define 
the entanglement verification under partial information 
on the performed measurements. Section IlIII starts with 
a discussion about different assumptions on the measure- 
ments. In addition we provide some further notation and 
background knowledge about the entanglement criterion 
that we employ for our purpose. Section FlII El finallv con- 
tains the above mentioned results for the two-qubit case, 
whereas Sec. IIVI is devoted to the general scenario of n 
uncharacterized dichotomic qudit measurements. Finally 
we conclude and comment on possible further extensions 
and directions in Sec. El Some technical details of the 
proofs can be found in appendix. 

II. PROBLEM DEFINITION 

Suppose that Alice and Bob observe an outcome prob- 
ability distribution denoted as Pr(a;, y|a, 6), where a la- 
bels different measurement choices with respective out- 
comes X for Alice and similar for Bob. These observed 
data have a quantum mechanical representation if there 
exists a quantum state pab and corresponding mea- 
surements, i.e., set of positive operator- valued measures 
(POVM) for Alice and {M^} for Bob such that 

Pr(a;, y|a, h) = tr(pABM^ ® M^), Vx, y, a, b. (1) 

The observed data are said to verify entanglement if and 
only if all states pab that satisfy this relation are en- 
tangled. Note that the measurement description is very 
crucial here because it ties a quantum mechanical mean- 
ing to the classical outcomes. 

In the following we consider the alternative that only 
partial knowledge is possessed on the measurement de- 
scription, meaning that the POVMs describing the mea- 
surement are not known completely. Hence each lo- 
cal measurement characterization is only assured to lie 
within a certain class. This set of possible POVMs will 
be denoted by A4a for Alice and Mb for Bob. In this 
case, successful entanglement detection implies that for 
all measurement descriptions only entangled states give 
rise to the observed data. More precisely, if S denotes 
the set of states having a quantum representation in ac- 
cordance with the assumed measurement description, 

S = {mb|3{M»} e Ma, {M'y} G Mb : (2) 

Pr{x, y\a, b) ^ tr{pABM^ ® Ml;), yx,y,a,b}, 

then the observed data Pr(2;, y\a, b) certify entanglement 
if and only if all state of this set S are entangled. 

Let us comment on the two extreme cases: If the mea- 
surements are completely specified each set M only con- 
sists of one possible element. In such case the question 
whether given observations verify entanglement is com- 
pletely answered for example with the help of entangle- 
ment witnesses, even if the measurement does not pro- 
vide full tomography (2^. In the other extreme that 



3 



the measurements are completely unspecified the sets M 
consist of all possible POVMs in all possible dimensions 
and the data exclusively correspond to entangled states 
if and only if a Bell inequality with the specified number 
of settings and outcomes is violated jl^ . 

Finally let us stress one more technical point: We do 
not assume any "convexification" of the problem as for 
example employed in Ref. ^28i | for dimension witnesses. 
Convexification would mean that if two different observa- 
tions Pi , P2 would have a separable quantum representa- 
tion then also its convex combination APi + (1 — A)P2 for 
all A g [0, 1]. However the problem is that the quantum 
representations might need different measurements, so 
that directly taking the convex combination on the level 
of quantum states does not work anymore. In Sec. IIII El 
wc provide an explicit counterexample, where convexifi- 
cation would lead to spurious entanglement detection. 



dimension restriction seems to us as the first non-trivial 
assumption that one can make [43j . 

Definition III. 1 (Qubit measurement models). For two 

dichotomic measurements, characterized by the operators 
A, A' according to Eq. (0) and satisfying Eq. we dis- 
tinguish the following cases: 

i) Qubit measurements: Both operators act on the 
same qubit. 

ii) Sharp qubit measurements: The POVM elements 
are rank-1 projectors on the same qubit, i.e., A, A' 
have eigenvalues ±1. 

Hi) Orthogonal qubit measurements: The eigenbasis of 
A and A' are mutually unbiased fJ2il- 



III. QUBIT CASE 



Parametrization of POVM elements 



This section concentrates on the two-qubit scenario. 
We start with the definition of different measurement as- 
sumptions followed by an explicit parametrization. Af- 
terwards we introduce the notation of a data matrix in or- 
der to express our results more compactly and also state 
the entanglement criterion that is employed to prove the 
main results in the last part of this section. 



In the following we introduce a parametrization of the 
POVM elements corresponding to different measurement 
scenarios. This parametrization will be convenient later 
for the technical proofs. Additionally it should further 
clarify the different measurement properties. 



Sharp but not orthogonal qubit m,easurem,ents 



A. Different measurement assumptions 

At first, let us specify the different measurement prop- 
erties more closely which were abstractly described by the 
set Ai in the previous section. We consider the simplest 
non-trivial case: Each party has two different measure- 
ment settings each of which having two different out- 
comes. Any of such dichotomic measurements is more 
compactly determined by the difference of two POVM 
elements, e.g., for the first setting of Alice we associate 
the operator 

A = AIIj^- AP^ ^ = i(l± A). (3) 



Here a; = ±1 labels the two different outcomes, while the 
resolution of the POVM elements M^^ follows because 
of normalization. In order that this operator A describes 
a valid POVM it must satisfy the conditions. 



yl-l>0, l-^>0. 



(4) 



Let us remark that this condition is still independent 
of any dimension restriction and it will reappear in the 
later section for the more general case of uncharacterised 
qudit measurements. The operator for the second choice 
of Alice is denoted by A' , while Bob's choices are given 
by B, B' respectively. 

The following definition summarizes the different qubit 
specifications that we consider. Let us point out that the 



In this case the operators A, A' can be written as fol- 
lows: 



A — cos(6')(Ti -f sin(6')o'j 
A' ~ cos(6')o'i — sin(6')o'j. 



(5) 
(6) 



where Ui , aj are two different possibly rotated Pauli oper- 



ators, i.e., they can be written as ai 



j-a and cr,- 



with two unit vectors 



S satisfying iii 



= 0. 



The parameter 9 characterizes the tilt between the mea- 
surement directions. Note that this relation can also 
be reversed, i.e., to express the Pauli operator in terms 
of the considered measurement operators. Formally the 
relation between orthogonal and non-orthogonal observ- 
ables is described by 



(7) 



" 1 " 




' 1 " 




= m 


A 






A' 



with kS 



R{0) 



2cos(6l) 
2cos(e) ~ 2sin(e) 



1 

2sin(e) 



(8) 



If we only refer to the 2x2 sub-matrix, formed by the 
second and third column and row, we use the label R2{0). 
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2. Non-sharp hut orthogonal qubit measurements 



D. Employed entanglement criteria 



When the measurements are not sharp but orthogonal, 
we directly employ the reverse parametrization 



(Ti = .Til + X2A 



(9) 
(10) 



with Xi £ R. Note that in order that A, A' correspond to 
physical observables given by Eq. ^ these parameters 
must satisfy X2 > 1 + and X4 > 1 + |x3|. Here we 
can choose without loss of generality X2,X4 > to be 
positive, by selecting appropriately <7i or — cr^. Formally 
sharp and non-sharp observables can be related by 



(11) 



" 1 " 




' 1 ' 






A 






A' 



with 



Six) 



1 

Xl X2 

X3 X4, 



(12) 



3. Uncharactensed qubit measurement 

The remaining case of totally uncharacterised qubit 
measurements can be considered as a combination of the 
above two cases. The overall transformation is given by 



1 



Rie)s{x) 



1 

A 
A' 



(13) 



The first operation S turns the operators A, A' into 
sharp, but not necessarily orthogonal measurements, 
which is considered afterwards by applying the transfor- 
mation R. 



C. Data matrix 

In order to express our results let us define some further 
notation. The observed data Pr(a;,y|a, 6) are compactly 
expressed in terms of a data matrix D3, given by the 
matrix of expectation values 



(1) (B) {B') 
(A) (AB) {AB') 
{A') (A'B) {A'B') 



(14) 



For convenience we often abbreviate the 2x2 sub-matrix 
containing only the full correlations, i.e., built up by the 
second and third row and column, as D2. Our criteria 
are typically given in terms of the singular values of this 
sub-matrix, denoted as X1/2 > 0. 



For entanglement detection we employ a criterion, 
which is a direct corollary of the computable cross-norm 
or realignment (CCNR) criterion [2^[33l. The corollary 
is formulated in terms of the singular values of the cor- 
relation matrix T3 given by 



(1) 



(afaf) (afaf) 



(15) 



Note that this correlation matrix T3 represents a special 
data matrix for which the employed measurement op- 
erators are sharp and orthogonal. Because of those sim- 
ilarities we employ the similar label T2 in order to refer 
to the full correlation 2x2 sub-matrix. 

Proposition III.l (Corollary of the CCNR criterion). 

Given the correlation matrix T3 with ordered singular val- 
ues Ao > Ai > A2 > 0. Then the CCNR criterion implies 
that any separable state necessarily satisfies 

IIT3II1 =Ao + Ai+A2 <2. (16) 

// the correlation matrix has vanishing marginals, i.e., 
(ct^) = {<^f) = for all k E {i,j}, and Aq = 1 then this 
condition is also sufficient. 

For completeness we provide a proof of this proposition 
in App. [X] With this stage set we will state in the next 
section our main results on entanglement verification in 
different two qubit scenarios. 



E. Main results 

In the following we state and prove our main results 
for qubits. We first consider the special case that the 
observed data matrix D3 has vanishing marginals. We 
obtain a complete solution for different scenarios if one 
only uses the knowledge of the appearing singular values, 
i.e., the criteria are minimized over all data matrices with 
fixed singular values. Using additional structure of the 
observation improves the detection strength as we will see 
later in Prop. IIII.31 For comparison reason the different 
detection regions are visualized in Fig. [T] 

Proposition III. 2. (Data matrix with zero marginals) 
Civen a data matrix D3 of the following form 



1 



Do 



(17) 



where the full correlation data matrix D2 is characterized 
by the singular values Ai > A2 > 0. These data verify en- 
tanglement under the assumption that both measurement 
of Alice and Boh are 

1 ) Sharp and orthogonal: Ai + A2 > 1 
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2) Sharp, non- orthogonal: + \/^ > 

3) Unsharp, orthogonal: Ai + A2 > 1 
4-) Qubit measurements: \/Xi + \/X2 > 



The definition of these properties is given in Def. \III. 1\ 
and these bounds are tight for the considered scenario. 

Remark III.l. Note that we assume that D3 actually 
originates from a quantum state under the considered 
measurement scenario, which can be assured for exam- 
ple if the singular values satisfy 1 > Ai > A2 > 0. 

Proof. Case 1) of sharp and orthogonal measurements 
has already been discussed in Ref. [1^. Alternatively it 
is a direct application of Prop. IIII.ll 

All other scenarios are proven along the following lines: 
Given the data matrix one first reconstructs the corre- 
sponding correlation matrix T3 by the appropriate trans- 
formations S{x),R{6) as given in Sec. IIIIBI In order to 
certify entanglement one employs Prop. IIII.ll However, 
since T3 depends on the transformation parameters, e.g.., 
9,x, . . . , one needs to optimize over all such choices. This 
will in general result in lower bounds on the singular val- 
ues of the data matrix. If these bounds are tight then the 
provided condition is necessary and sufficient in order to 
detect entanglement with provided data. 

Case 2) At first let us concentrate on the sharp but 
non-orthogonal case, in which the correlation matrix is 
given by T3 = R{a)D-iR{f3y . For the block-diagonal 
data matrix the resulting correlation matrix is of similar 
block structure, i.e., T3 = diag[l,T2] with 



T2 = R2{a)D2Rl{P). 



(18) 



Hence if the ordered singular values of T2 are denoted 
as ti > t2 > Q then Prop IIII.ll states that the state is 
entangled if and only if ti -I- ^2 > 1 holds for all values 
of a, /3. In order to minimize the sum ti + t2 over the 
angles we first lower bound this quantity by an expres- 
sion containing only the singular values of the appearing 
transformations since this is more easily optimized in the 
end. 

The lower bound is derived using the inverse relation 
of Eq. (UHl), 



D2 = R2{a)-^T2Rl{P)- 



with 



R2{a)- 



cos(a) sin(a) 
cos(a) — sin(Q;) 



(19) 



(20) 



In the following discussion we employ the abbreviations 
oi > 02 > and 61 > 62 > for the ordered singular 
values of R^{q)~^ and R'^iP)''^ respectively. Further- 
more recall that D2 is characterized by its two singular 
values Ai > A2 > 0. For the matrices on the left- and 



right-hand side of Eq. (IT9|) the following relations hold 

Ai A2 



tl + tl 



> 



aia2bib2 ' 

(Al+A2)2 



afbl 



ajbl 



(21) 
(22) 



The proof of these two relations involves some technical 
details and is given in App. IB II Employing these two 
identities provide 



[ti + t2f > min 

a,/3 



(Ai-f A2)2 



A1A2 



afbl 



albl 



01026162 



> 



Ai + VA^ 



(23) 
(24) 



The last inequality arises if one employs the true singu- 
lar values oi (a) , . . . and performs the minimization; for 
an explicit proof of this optimization see Lemma [T] in 
App. O Eq. p4|) confirms that the state is entangled 
if and only if -y/AT -I- \/Xi > v^, where the sufficiency 
follows from the fact that all appearing inequalities can 
also be achieved with equality. This finilizes the proof of 
the Case 2). 

Case 3) Next consider the unsharp, orthogonal case. 
The correlation matrix T3 = S{x)D3S{y)'^ , with ap- 
propriate "sharpener" transformations S{x) given by 
Eq. (fT2)) . has the following block structure 



T3 = 



1 




■ 1 




'1 y^' 


X Sx 




^2 







X xy'^ + SxD2Sl 



(25) 
(26) 



Here we used an analogue block decomposition for S{x) 
with column vector x — [xi^x^Y' and diagonal sub- 
matrix Sx = diag[2;2, X4\ and respective abbreviations for 
the transformation of Bob. In order to prove the state- 
ment we will use of the following three inequalities for 
the ordered singular values of , 



h > 1, 
totit2 > A1A2. 



(27) 
(28) 
(29) 



The proof of this inequalities is given in App. IB 21 Fur- 
ther, as shown in Lemma [2] of App. [C] these conditions, 
together with the ordering condition 1 > Ai > A2 > 0, 
ensure 



tQ+ti+t2> I + A1+A2. 



(30) 



If this ordering is not valid, i.e., Ai > 1, entanglement 
directly follows because of 



to + ti+t2>to + ti> > 2 



(31) 



via the inequality of arithmetic and geometric means. 
Thus in total Ai -f- A2 > 1 is necessary and sufficient for 
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entanglement, sufficiency because one detects the same 
as in the more restrictive case of sharp, orthogonal mea- 
surements. This finishes the proof for the non-sharp, 
orthogonal case of qubit measurements. 

Case 4) For the remaining scenario of fully uncharac- 
terised qubit measurements we can largely employ the 
previous results. In this scenario the correlation matrix 
T3 = R{a)S{x)DsSiyfR{(3 f is given by 



1 



i?2(a) 



1 

X xy'^ + SxD2S, 

~T 

y 



X xy'^ + R2{a)S,D2S^R2if3f 



(32) 
(33) 



with X = R2{a)x, y = i?2(/3)y. In this case the important 
sub-matrix is 

T2 = R2{a)S^D2S^R2Wf = R2{a)D2R2Wf, (34) 

which can be considered as the central sub-matrix of the 
sharp but non-orthogonal case, Eq. (fT8|) . but where the 
transformation is applied to D2 instead of the true data 
matrix D2 itself. From the sharp, non-orthogonal case we 
know that the singular values of this matrix T2 , denoted 
as ti > t2 satisfy 



tl+t2> 



1 



Ai + 



2 1 



A1 + VA2 , (35) 



where Xi are singular values of 1)2. Next, using similar 
arguments as already presented in the sharp but non- 
orthogonal case we can derive the following set of in- 
equalities for the singular values of T3: 

to > 1, (36) 
toh > fi, (37) 
totit2 > hh- (38) 

Using once more Lemma [2] and Eq. ([35|) provides 

1 / / , \2 

to+ti+t2>l + ti+t2>l + - V Ai + VA2 j (39) 

if one has the ordering 1 > h > t2 > 0. If > 1 one ver- 
ifies entanglement again by the inequality of arithmetic 
and geometric means. This finally shows that the state 
is entangled if and only if + y/X2 > which proves 
the claim for uncharacterised qubit measurements. □ 

Next let us provide an important numerical example. 
First it demonstrates that the unsharp, orthogonal case is 
indeed different from the completely characterized case. 
Together with the previous proposition it shows that the 
sharp, non-orthogonal and unsharp, orthogonal case are 
indeed inequivalent to each other, i.e., there are observa- 
tions which are exclusively detected by one of these two 
scenarios. Additionally, this example proves that the en- 
tanglement verification in the unsharp case is in fact a 
non-convex problem. This means that one must be very 
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FIG. 1: DifTerent detection regions for a data matrix with 
singular values Ai, A2. The solid and long-dashed lines corre- 
spond to the case of a data matrix with vanishing marginals 
as discussed in Prop. IIII.2I The remaining two lines corre- 
spond to determinant detection rule det(_D) — A1A2 given in 
Prop. ITVTTI for qubits (dotted) and qutrits (short-dashed). 



careful in applying Prop. IIII.21 it is for example not pos- 
sible to use it on a "depolarized" version of the observed 
data matrix D3, i.e., the one that one obtains by setting 
the marginals equal to zero. 

Example III.l. The data matrix 
1 1 - V3 



D. = 



l-VS (15-8V3)/2 



1/2 



1 

-0.73 



-0.73 
0.57 



0.5 



(40) 



(41) 



can originate from a separable state in the case of un- 
sharp orthogonal measurements, but verifies entangle- 
ment for .sharp, non- orthogonal measurements. 

Moreover it shows that the unsharp scenarios are in- 
deed non-convex problems, i.e., convex combination of 
two separable data matrices might not be separable any- 
more. 

Proof. First let us give the separable state and its cor- 
responding measurements that are consistent with the 
given data matrix. This also discloses the generation of 
this example. The data matrix given by Eq. (|40)) is ob- 
tained by measuring the separable state 



Pscp 



1 (g) 1 -I- i ((Tj; ® cr^ + fj^ ® 0-2 



(42) 



with unsharp measurements, A parametrized according 
to Eq. (jni) with ai = a^, Xi = 1 -\- y/3 and X2 = 1 -\- xi, 
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while A' = (Jz and the same settings for Bob's side. The 
reason for the choice of Xi becomes clear afterwards. 

Verifying entanglement under the premise of sharp 
qubit observables goes as follows: Note that if the mea- 
surements are even orthogonal applying the CCNR cri- 
terion of Prop. nil. II immediately shows entanglement. 
In the non-orthogonal case one utilizes again the trans- 
formations R{a),R{p) in order to generate T3 as done 
in the previous proposition. Even with non-vanishing 
marginals the important sub-matrix is given by T2 — 
i?2(o)^2-R2(/3)^ with Z?2 being the sub-matrix contain- 
ing the full correlations. Though not directly stated as 
the CCNR criterion, a state is already entangled if the 
singular values of T2 satisfy ti+t2 > 1 H^. As shown in 
the proof of the previous proposition this relation is as- 
sured if the singular values of D2 fulfil y/Xi + ^/X2 > V^- 
Since the data matrix given by Eq. (|40|) satisfies this con- 
dition this proves that cannot be compatible with a 
separable state under sharp qubit measurements. 

Finally one needs to verify that the data matrix given 
by Eq. (|40p is at all consistent with a valid quantum state 
using sharp measurements. This is necessary in order to 
assure that the set S defined as in Eq. ([2]) is indeed non- 
empty. However the operator 



Pont 



ycTy ®ay+ ^ [D■i\^jai ® ajj (43) 

i,je{0,x,y} 



with y = 4\/3 — 7 represents a valid state compatible with 
the data matrix if one employs the sharp measurements 
A = B = ax, A' = B' = a^. Let us point out that this 
physicality condition determined the parameter xi: We 
optimized the detection condition y/Xi + while still 
keeping the data compatible with a valid state. 

Furthermore, this example provides an explicit in- 
stance for the failure of convexity [l^. If D3 corre- 
sponds to a separable state then also its marginal in- 
verted version because it effectively only represents 
a classical outcome interchange -1-1 -H- —1 on both 
sides. Thus also the data matrix given by Eq. PO]) with 
marginals —(1 — VS) is separable. However, taking the 
equal mixture leads to the "depolarized" data matrix 
D3, = diag[l,L'2] which would verify entanglement ac- 
cording to Prop. inL2] □ 

The following proposition demonstrates that the en- 
tanglement properties are not fully determined by the 
singular values of the observed data matrix. Moreover, 
for this special data structure it is interesting to observe 
that the extra knowledge of sharpness and orthogonality 
is irrelevant for the detection strength and mere informa- 
tion about the dimension of the measurements suffices 
to verify the same fraction of entanglement. Further- 
more, since these observations satisfy all CHSH inequal- 
ities [10[, complete device independent detection is not 
possible. 

Proposition III. 3. (Diagonal data matrix) Observa- 



tions of a diagonal data matrix 

1 



Ai 



A2 



(44) 



with 1 > Ai/2 > verify entanglement under the assump- 
tion of qubit measurements if and only if Xi + X2 > 1. 
Hence one verifies the same fraction as with sharp, or- 
thogonal qubit measurements. In contrast, complete de- 
vice independent entanglement verification fails. 

Proof. The proof runs analogous to the qubit measure- 
ment scenario of Prop. IIII.21 Note that the full cor- 
relation matrix T3 is given by Eq. (j33[) and that one 
detects entanglement if and only if the singular val- 
ues of T2 = R2{a)D2R2{PY given by E q. ((M| . satisfy 
ii + ^2 > 1- However in contrast to Prop. IIIL2] it is now 
possible to derive a tighter lower bound by exploiting the 
diagonal structure of 



D2 



X2y2Xl 




■ Ai 


y2yi^2 




A2 



(45) 



The singular values of D2 are given by the diagonal en- 
tries, which fulfil Ai > Ai, due to the constraints on the 
parameters Xi,yi. In order to finish the proof we em- 
ploy as usual certain inequalities for singular values of 
the matrices T2 and D2: 



tit2 > A1A2, 



(46) 
(47) 



These inequalities imply ti -\- 12 > Xi + X2 and therefore 
prove the claim of the proposition. 

In order to prove Ineq. PS|) one applies the deter- 
minant multiplication rule together with the property 
I det[i?2(')]l ^ 1 th^-t '^^^ be checked directly from the 
definition given by Eq. ([8]). The second ineq. (|47| is ver- 
ified by 



P1+P2 



tr(T2Tj) 



(48) 



= — |(Ai + A2)2[csc(Q)2csc(/3)2H-sec(a)2sec(^)21 
16 



+ (Ai-A2)^[csc(a)2sec(/3)2-t-sec(a)2csc(/3)2] } (49) 

>l [(Ai + + (Ai - A2)2] >Xl + Xl, (50) 

where the first inequality is obtained by minimizing each 
term within the squared brackets separately. 

Complete device independent entanglement verifica- 
tion fails because all possible Bell inequalities are sat- 
isfied, which is equivalent to a separable quantum repre- 
sentation in the bipartite case [iB,|il. □ 



IV. QUTRITS AND BEYOND 

As shown in the previous section solely having the 
knowledge that one is measuring a qubit is enough to 
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detect entanglement even if the corresponding Bell in- 
equalities, the set of inequivalent CHSH inequalities, are 
not violated. Nevertheless if all these Bell inequalities 
are satisfied then the observed data can be reproduced 
by appropriate measurements onto a higher dimensional 
separable state for the considered case this would 

be in dimensions 4 (g) 4. Thus the only other non-trivial 
case is the instance of qutrits. In the following we prove 
that even the qutrit assumption alone suffices to detect 
more than with the CHSH Bell inequalities. 

Rather than defining different notions of sharpness 
or orthogonality for higher dimensional measurements 
or more settings, we focus on the completely uncharac- 
terised case of n dichotomic measurements on a c? dimen- 
sional system. Each dichotomic measurement is uniquely 
determined by the operator given by the difference of two 
POVM elements and is denoted as Ai with i = 1, . . . ,n 
in the following. The only defining inequality for all 
these operators Ai, besides that they are all acting on 
the same c?-dimensional Hilbert space (D'', is the condi- 
tion of Eq. (U) which ensures that they correspond to 
valid quantum measurements. Similar conditions are im- 
posed for the measurements for Bob labelled as Bi. We 
employ once more the notion of a data matrix D, each 
entry defined as [D]ij = {Ai (X) Bj) for i,j = 0,...,n, 
with Aq = Bq = 1, such that it also contains the ob- 
served marginals. Obviously we also have [D]oo = 1, 
which we always assume to be fulfilled if we speak about 
a data matrix. This is again employed to provide a more 
compact solution, which is stated in the following propo- 
sition. Besides from the above mentioned qutrit example 
in the CHSH case it has a few more consequences which 
are commented afterwards. The condition for qubit and 
qutrits are plotted in Fig. [T] 

Proposition IV. 1. (Data matrix for n dichotomic mea- 
surements on two qudits) If the data matrix D corre- 
sponding to n dichotomic measurements satisfies 



\det{D)\ > 



1 



n+l 



(51) 



one verifies entanglement under the assumption of d- 
dimensional measurements. 

In case one has at least as much settings as dimen- 
sions, i.e., n > d, already the condition 



|det(L>)| > 



d-1 



(52) 



ensures entanglement. For the special case of two settings 
n = 2 and qutrits d — 3 the condition can be improved to 



|det(D)| > 



64 
81' 



(53) 



Proof. We follow a similar proof technique as in the pre- 
vious section. The data matrix is transformed with ap- 
propriate corrections Ga,Gb in order to form a kind of 



correlation matrix C = GaDG^ for which one employs a 
known entanglement criterion. 

This correlation matrix C is very similar to the previ- 
ously employed matrix T^. Here it is defined as C[p]ij = 
tr[pKf- ® K^) where each local set Ki consists of or- 
thonormal (with respect to the Hilbert Schmidt inner 



product) observablcs, i.e., iv{KiKj) 



Using the in- 



clusion principle in a similar fashion as in the proof of 
Prop. IHI.ll one finds that the state is entangled if the 
singular values of C, denoted as Ci, fulfil 



(54) 



In the following we bound this trace norm by the de- 
terminant of the correlation matrix | det(C)| = Ci and 
the extra knowledge that the largest singular value satis- 
fies Co > 1/d which is implied by the choice Kq ~ 1/ \fd 
and the inclusion principle. This provides the following 
estimate, 



> Ci > min Cn + > c^ > min Cn + n | [ C; (55) 

^ CO>l/d ^ co>l/d \^^. I 



> min Co 

co>l/d 



co>l/d 

det(C)|\ " 



Co 



'(n + l)|det(C)|^ if|det(C)|^ > i 
i +n(d|dct(C)|)" else 



(56) 



(57) 



Here we employed the inequality of arithmetic and geo- 
metric means in the second step, while the optimization 
given by Eq. (j56p is performed using standard analysis. 
Note that the first solution in Eq. ([57)) is the uncon- 
strained optimum that could have been inferred directly 
by applying the inequality of arithmetic and geometric 
mean to all terms. However, under the constraint on the 
largest singular value this solution is only reached if the 
determinant of the correlation matrix satisfies the stated 
extra condition. This distinction is necessary for the im- 
proved condition in case n > d. 

The remaining strategy is to lower bound each solu- 
tion of Eq. ([57)) by an expression involving the data ma- 
trix. Afterwards one investigates which conditions assure 
that this lower bound actually exceeds one such that, in 
spirit of Eq. ([54)) . it would signal entanglement. The 
more stringent of these two cases will be the final entan- 
glement criterion. Here we need the following inequality 
that relates the determinant of correlation and data ma- 
trix. 



|det(C)| 



> 



|det(D)| |det(Ga) 
|det(i^)|d-("+i) 



\det{Gl) 



(58) 
(59) 



which follows from the bound | det(G)| > d 2 for each 
of the above mentioned transformation Ga, Gb, proven in 
App.jDj 



9 



Let us start with the second solution and employ 
Eq. ([5^ which results in 



3+n(d|det(C)|)" > 3 (l + ?^|detp)|- 
a a 



(60) 



which is larger than one if and only if the condition given 
by Eq. ([5^ holds. For the first solution one obtains 



(n + f)|det(C)|— > 



1 



1, |det(i:»)|" 



(61) 



The first part in the maximum follows from the region 
constraint on | det(C)|, while the second part is obtained 
using Eq. (|59p. The maximum appears because both 
bounds are valid. The right hand side of Eq. (|6T|) is larger 
than one if already one the terms does so. For the general 
case one chooses the second part of this maximum, which 
is larger than one if and only if the determinant of the 
data matrix satisfies Eq. ([51]) . Because this condition is 
weaker than the previous condition from Eq. this is 
the entanglement criterion for the general case. For the 
special configuration of n > d we employ the first part 
of the maximum since it always exceeds one. Hence only 
the condition from Eq. (|60p is the relevant for this case, 
which proves the part of the special case n > d. 

The improved condition for the qutrit case follows 
from a sharper lower bound on the transformation Ga , Gb 
given by |det(G)| > ^3/8 which is proven in App. |DJ □ 



It is worth to stress that one can employ Prop. ITvTll not 
only with the determinant of the full data matrix D, but 
also for each sub-determinant. This describes the case 
that certain measurement settings are left out, which is 
useful when two or more settings coincide or are linearly 
dependent in which case the determinant of the whole 
data matrix vanishes. 

Interestingly Prop. IIV.II detects bound entanglement. 
In what follows we provide an explicit example of a PPT 
bound entangled state that is detected via the criterion 
given by Prop. IIV.II i.e., solely by having knowledge 
about the underlying dimension. The state 



PBFP = ^(*iB*A': 



*to*t,o, +1'7o3'7,o, (62) 



AB ^ A'B' 



AB^A'B' 



+ *AB*A'B' + *AB*A'B' + *AB*A'B') 

with $+, . . . , denoting the projectors onto the stan- 
dard two-qubit Bell states, has been shown to be 4 ® 4 
bipartite PPT bound entangled [3l| under the splitting 
AA'jBB'. Assume that one performs "good" enough 
measurements described by all traceless operators Aki = 
® f/^ built up by tensor products of the identity and 
the Pauli operators, and the same measurements for Bob. 
Mixed with white noise p{p) = (1 — p)pbfp +pl/16, the 
criterion given by Eq. ([52]) becomes 



\dem\ = 



15 



> 



15 



(63) 



and thus verifies entanglement as long as p< 2/5 = 0.4. 
This values seems to coincide with the point where en- 
tanglement disappears, i.e., for p > 0.41 the state is sepa- 
rable using the method of Ref. [33l ■ This detection capa- 
bility represents a clear advantage over Bell inequalities, 
where it is still unknown if measurements on a bipartite 
bound entangled state could at all violate a Bell inequal- 
ity, though recent results seem to falsify this belief [33| . 

To conclude this section we consider the two-qubit 
Werner states pw{p) = (1 — p)"^" +pl/4 and three or- 
thogonal standard measurements, i.e., ax,cFy, cTz for both 
sides. Then the above stated conditions verify entan- 
glement as long as the white noise parameter satisfies 
p < 2/3 ~ 0.67, which coincides with the point where 
entanglement vanishes. This represents another impor- 
tant improvement with respect to Bell inequalities, since 
first there is hardly any good Bell inequality known that 
detects entanglement if p > 1 — l/-\/2 « 0.29 and second, 
above p> 7/12 w 0.58 it is known that no measurement 
would violate a Bell inequality [sj]. 



V. CONCLUSION & OUTLOOK 

We have investigated the task of entanglement detec- 
tion for cases where only some partial information about 
the performed measurements is known or assumed. The 
considered scenarios included properties like sharpness, 
orthogonality or where only the dimension of the underly- 
ing measurements is fixed. Via this extra information one 
verifies more data as resulting only from entangled states 
than in the totally device independent setting while still 
keeping a good detection strength in comparison to the 
fully characterized case. 

There are many further research lines connecting from 
here: A thorough investigation of higher dimensional 
states and more measurement settings is clearly interest- 
ing in order to clarify the power but also the limitations 
of this intermediate approach. For that one should be 
aware that the current methods only represent first tech- 
niques towards these directions. Here alternative tools 
might be necessary, even an efficient numerical approach 
would be of great help. Another approach would be the 
investigation of detection methods for the multipartite 
case in a similar intermediate setting. Since our criterion 
rests on an entanglement criterion based on the corre- 
lation matrix this extension might be possible utilizing 
recent detection methods for genuine multipartite entan- 
glement using the correlation tensor [35l [36j . Regard- 
ing explicit tasks, since our results show that one veri- 
fies even entanglement if the underlying Bell inequality 
is not violated this means that the lower bound on the 
concurrence given in Ref. (37j could be improved. This 
partially characterized scenario might also be useful in 
order to obtain steering equalities which are more robust 
against calibration errors, in similar spirit as in Ref. [ssj . 
Finally it is tempting to apply these result also to quan- 
tum key distribution operated in a similar intermediate 
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setting as described here, which has already been started 
in Ref. [s^. For that in particular Prop. lIlL3l is interest- 
ing because it describes exactly the kind of observations 
that one expects in an entanglement based BB84 proto- 
col. Since it states that one verifies the same fraction 
of entanglement as with totally characterized measure- 
ments, this hints that a very strong "semi-device inde- 
pendent" key rate could be obtained if one just possesses 
the knowledge that one measures a qubit. This would al- 
low a much larger freedom in finding appropriate squash 
models for quantum key distribution since the measure- 
ments don't have to be fixed anymore [2]| . 
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Appendix A: Proof of Prop. IIII.ll 

In this appendix we provide the proof of the CCNR 
criterion adapted to our partial information setting. 

Proof. It is only necessary to consider the case of partial 
information, since the CCNR criterion is typically for- 
mulated in terms of the full density operator. Following 
Ref. the CCNR criterion can be expressed as follows: 
Let T4 denote the complete correlation matrix of two 
qubits, which results from T3 by adding the remain- 
ing Pauli operator for each local side. Suppose that 
its corresponding ordered singular values are denoted as 
to > h > t2 > ts > 0. Then the CCNR criterion states 
that for any separable state these singular values fulfil 

According to the inclusion principle, cf. Corollary 3.1.3 
of Ref. [4l| , the ordered singular values are lower bounded 
by the singular values of any sub-matrix. Thus one ob- 
tains, ti > Xi for i = 0, 1, 2, such that one arrives at 

3 

Ao + Ai + A2 < < 2. (Al) 

i=Q 

Whenever this condition is violated the state must nec- 
essarily be entangled. 

In case of vanishing marginals with Ao = 1 the con- 
dition transforms into Ai -I- A2 < 1. Note that the pa- 
rameters Ai/2 < 1 are also the singular values of the 
sub-matrix T2. Using similar techniques as presented in 
Ref. [13 it is possible to fit a rotated Bell-diagonal sepa- 
rable state to these data. This is achieved along the fol- 
lowing lines: First, consider another correlation matrix 



T[p\ai3 = tr(p(TQ (g) cr^) which is built up by the standard 
Pauli operators {ax^Cy^az} for each local side. Assume 
that the given sub-matrix T2 is precisely the upper-two 
block of such a correlation matrix, i.e., T = diag[r2jO] 
filled with additional zero entries. The singular value de- 
composition of this correlation matrix is given by T = 
OaAOb with singular values A = diag[Ai, A2, 0] > and 
Oa, Ob being special orthogonal matrices of similar block- 
diagonal form, i.e., Oa = diag[Oa,±l] and an analogous 
form of Ob. Here note that Oq, Ob are the, not necessar- 
ily special orthogonal matrices form the singular value 
decomposition of T2 = Oa diag[Ai, A2]0^. 

Next let us discuss the special case of a diagonal cor- 
relation matrix, i.e., T[p\ = diag[Ai, A2, 0]: These data 
correspond to a Bell-diagonal state, abstractly expressed 
as pbs = X^jPi {hsi\ with standard Bell states \bsi) 
and appropriate weights equal to = (1 ± Ai ± A2)/4 
having all four combinations. The above stated condition 
ensures that all probabilities are indeed non-negative and 
are upper bounded by 1/2, which ensure separability in 
this case [42| . 

For the general case one employs the relation that any 
special orthogonal transformed correlation matrix cor- 
responds to some special unitary transformation on the 
level of quantum states [i^ , 

OaT[p]Ol - T[Ua ® UbpUl ® Ul] . (A2) 

Since local unitary transformation do not change the 
entanglement properties, the appropriately transformed 
state Ua ® UbphsUl ® Ul is the actual separable state for 
the general correlation matrix T. Finally, this rotation 
property is exploited once more to support the assump- 
tion that T2 is the sub-matrix of the first two rows of T , 
since any two orthogonal vectors can be rotated such that 
it matches these axis. This finally proves the claim. □ 

Appendix B: Details for two qubit case 

1. Relations for sharp nonorthogonal 
measurements, Eqns. (|2ip . (|22|) 

Equation (PT|) is a direct consequence of the determi- 
nant multiplication rule, i.e., dei{AB) = det(yl) det(i?), 
and that the absolute value of the determinant is equal 
to the product of its singular values. In order to de- 
rive the second inequality one employs the singular value 
identities 

fc fc 

Y,<^^iAB)<Y,<^^iA)<J^{B), (Bl) 

i=l 1=1 
fc k 

Y[a,{AB)<Y[a, {A)a, {B) , (B2) 

1=1 i=l 

where (Ji(-) denotes the decreasing ordered singular val- 
ues, cf. Ref. [4l|. Applying these identities to Eg. [TOl 
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leads to 

Ai + A2 <hiai{[R2{a)]-^T2) + &2CT2([-R2(a)]"'r2) (B3) 

<6l(7l([i?2(a)]"'T2) 

+ h2[aiti + a2t2 - (Ti(i?2(a)"^T2)] (B4) 
<{aibi)ti + (a2^>2)t2 (B5) 

<^{alhl + albl){tl+q), (B6) 

where the Cauchy-Schwarz inequality is applied in the 
last step. 



Relations for non-sharp orthogonal 
measurements, Eqns. (|27[l - (|29|) 



Lemma 1. 



The first conditiongiven by Eq. \21\ follows from the 
inclusion principle j4ll | using the first entry as a sub- 
matrix. 

The last inequality Eq. (PH]) holds because of the de- 
terminant multiplication rule, 



toht2 = I det[5(x)]|| detp3)|| det[5(y)]| 



(B7) 



= I dct{S.,D2Sl)\ = (.T2a;4)(AiA2)(2/2y4) (B8) 
> A1A2, (B9) 

and the bounds on the appearing parameters, e.g., X2 > 

1 + I.Tll > 1. 

In order to prove Eq. psp we apply the inclusion prin- 
ciple to a particular chosen 2x2 sub-matrix. For this ar- 
gument the matrix D2 = SxD2Sy attains special impor- 
tance. First, note that the singular values of D2, denoted 
as Ai > A2, satisfy A^ > A^ because the transformations 
satisfy Sx,Sy — 1 > 0. Next, employ the singular value 
decomposition D2 = U'EV'^ with S = diag[Ai, A2]. Since 
the singular values of T3 remain invariant under orthog- 
onal transformations, we apply appropriate orthogonal 
matrices to diagonalize D2, which leads to 



1 



1 



V 



1 

-T - -T 

X xy - 



(BIO) 



with X = U"^x,y = V^y. Using the sub-matrix formed 
by the first two rows and columns one obtains 



1 2/1 _ 

xi xiyi + Ai 



(Bll) 



which has determinant Ai > Ai. Then the inclusion prin- 
ciple directly states Eq. 



Appendix C: Optimization problems 

In this appendix we prove two lemmas concerning opti- 
mization problems appearing in the proof of Prop. IIII.2I 
They are mainly given because of completeness of the 
manuscript. 



Suppose Ai, 
(Ai4 



A2 > 0. 

A2)' , 



Then the solution of 
r> -^1-^2 



aibl 



01026162 



(CI) 



with Oi = cos(a)2 > 02 = y^2sm{a)^ and similar for 
hi with another angle (3 is 



Ai + 



(C2) 



Proof. First note that the ordering of the singular val- 
ues does not modify the solution if the optimization is 
performed over the full period of each angle since wrong 
ordering only leads to larger function values. Using the 
parametrization a = ^ + 5,l3 — simplifies the func- 

tion to 



(Ai + A2)^ 



4A1A2 



2 -I- cos(47) -I- cos(4(5) | cos(47) - cos(4(S)| ' 



(C3) 



Hence it effectively only depends on u = cos(47) and 
V = cos(4(5) which are both in the interval [—1, 1]. Note 
that the boundary of this feasible set is characterized by 
either u, v taking on the extreme values of this interval. 
In the following we prove that the minima lies at this 
boundary. 

Using the linear variable transformation x = u + v, 
y ~ u ~ V changes the function to 



(Ai+A2)2 4A1A2 



\y\ 



(C4) 



Now, taking partial derivatives, one finds that this func- 
tion is decreasing with respect to x in the interval (—2, 2] 
and depending on the sign of y also decreasing in the +y 
OT —y direction [the exceptional cases y = ov x = —2 
can be excluded since they are no minima]. Since the 
variables x, y are bounded this shows that the minima 
is attained at the boundary. Going back to the form 
given by Eq. IC3I means that either cos(47) = ±1 or 
cos (4(5) = ±1. Here it does not matter which cosine 
is put to its extreme values such that one only needs to 
consider one of them. Only one of these values ±1 leads 
to the solution, for the other one can directly verify that 
its solution is larger than given by Eq. IC2I Concluding, 
only the following function needs to be optimized over 
the angle ■0, 



4A1A2 



(Ai + A2)^ ^ 

3 4- cos(V') 1 — cos(-0) 



(C5) 



This can be performed directly by looking for the van- 
ishing derivatives and using only its real solutions. This 
finally leads to the solution given by Eq. IC2I □ 

Lemma 2. Suppose Aq > Ai > A2 > 0. Then the solu- 
tion of 



mm Mo + A^i + M2 

s. t. Mo > Aq, MoMi > AoAi, 
MoMiAi2 > A0A1A2, 
Mo > Ml > M2 > 0, 



(C6) 
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is Ao + Ai + A2 



Appendix D: Transformation determinants 



Proof. In order to prove the lemma let us first state the 
solution of the following sub-problem, 



A 



mm X ~\ 

a;>a;min>0 X 



2VX 



if VX > X-n 

else 



(C7) 



with A > that follows using standard analysis. Note 
that \/A is the argument of the unconstrained minima. 

Now we turn to the intended problem given by Eq. IC6I 
At first consider the case that the parameter ^0 is fixed 
and that we bound the sum of the remaining two param- 
eters from below, which provides 



min IJ.1 + Hz > niin fii 

Ml,Ai2 fJ.l>fJ.l 



AoAiA2\ 1 
Ml 



> 



MO 



if ^J'0 > P-o 



^ + A2 if Ao < A^o < flo 



(C8) 
(C9) 



using the abbreviations p,i = AqAiZ/xq and p,Q = A0A1/A2. 
In the first inequality we employ the lower bound on ^2 
from the problem formulation. The second inequality is 
an application of the sub-problem in which the conditions 
are re-expressed in terms of the parameter fxo- In the 
following we use these lower bounds and consider vari- 
ations of /io within the corresponding valid region. For 
simplicity let us assume strict inequality Aq > Ai > A2; 
the statement with equality of some or all parameters 
follows then by continuity. 

Let us start with the case Aq < /-to < Ao- Using the 
derived lower bound gives 

min uq + ui + 112 > min uq H ^^—^ + A2 (CIO) 

AoS/ioSMo 



> Ao + Ai + A2, 



(Cll) 



via another application of the given sub-problem. Note 
that the unconstrained minima is not achieved, i.e., 
VAqAi Ao because of the strict ordering. 

Next consider the case iJ-o > p-o, for which we have to 
employ the other lower bound in Eq. IC9I This leads to 
the 



min uo + Ail + M2 > min uo + 2^ /_2^L^ (C12) 
Mi:M2_ ^" ' Mo> An V Mo 

Mo > Mo ' 

>Mi + 2A2>Ao + Ai+A2. (C13) 
A2 



The optimization problem appearing in Eq. IC12I is very 
similar to our sub-problem, in particular it is convex 
again for positive ^q. Its unconstrained minima is at 
the argument A0A1A2 which, however, is outside the 
allowed region, i.e., /2o > v^AoAiAi^ due to the strict or- 
dering of the A's. Thus the optima is attained at the 
boundary /io = Mo- The last inequality represents an- 
other consequence of the ordering property since it is 
equivalent to Ao(Ai — A2) > A2(Ai — A2). □ 



In this appendix we provide a proof for the bounds 
on the determinant of the transformations G used in the 
proof of Prop. IIV.II It is a direct consequence of the 
Gram-Schmidt procedure. 

Lemma 3. The linear transformation G that maps the 
operators from the data matrix {Ai} all acting on (C* 
with Kq = 1 and — 1 < Aj < 1 for all i = 1, . . . ,n into 
an orthonormal operator set {-ft',}, i.e., tT{KiKj) — Sij 
fulfils 

|det(G)| > d-"^. (Dl) 
For the case n = 2, d = 3 the bound can be improved to 



|det(G)| > 



V3 
8 ' 



(D2) 



Proof. The linear operation G can be obtained for in- 
stance using the Gram-Schmidt process. Without loss 
of generality this procedure can be decomposed into 
two operations G — ON. The first transformation N 
should map each operator Ai to its normalized form 
Ai = Ai/ ^ytI{Af), such that the second transformation 
O only needs to orthogonalize them. This second linear 
operation always stretches the "vectors" Ai such that the 
volume spanned by this set always increases, therefore 
|det(0)| > 1. The first transformation is a diago- 
nal matrix with entries given by 1/ ■\/tr(A|). Since each 
operator Ai is bounded by the identity due to the restric- 
tion that it describes a valid measurement, one obtains 
tiA^ < trl < fi or finally 

I det(G)| = I det(0)| I det(iV)| > ^"("+1)/^ (D3) 

For the special case of ?i = 2 and c? = 3 we explic- 
itly carry out the Gram-Schmidt process and minimize 
the determinant under the given constraints. We con- 
sider the case that G maps the operators to following 
orthonormal set {A'o = l/\/3, A'l, A'2}. This resulting 
operation G is of triangular form, i.e., 



G„ — 



1 

V3 

* A^i 

* * A2 



(D4) 



which has a determinant N1N2/ V3. In order to obtain a 
lower bound one needs to minimize N1N2 or respectively 
maximize l/(A'^iA2)^, that results in 



tv{A') 



MA)? 



triA") 



,2^ MAT 



< 



< 



tiiAA') 

tr{A^ 
2 



tr{A) tr(A') 



MA)] 



tT{A 



I2\ 



MA')] 



(D5) 
(D6) 
(D7) 
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The last inequality originates from the bound 
tr(A2 



,2^ MA)r 



- {Ij + q + ll - hh - hh - hh) < 3 (D8) 



constraint < 1 due to the measurement condition 
of Eq. m In this expression at most two of the last 
three terms can be positive but one must necessarily 
be negative. Suppose that this term is —hh < 0, 
then I2 + I3 — hh < 1 holds given the stated box con- 
straints. □ 



where k are the eigenvalues of A that satisfy a box 
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